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We investigate the effects of phase noise and particle loss on the dynamics of a Bose-Einstein 
condensate in an optical lattice. Starting from the many-body master equation, we discuss the ap- 
plicability of generalized mean-field approximations in the presence of dissipation as well as methods 
to simulate quantum effects beyond mean-field by including higher-order correlation functions. It 
is shown that localized particle dissipation leads to surprising dynamics, as it can suppress decay 
and restore the coherence of a Bose-Einstein condensate. These effects can be applied to engineer 
coherent structures such as stable discrete breathers and dark solitons. 
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I. INTRODUCTION 

Decoherence and dissipation, caused by the irreversible 
coupling of a quantum system to its environment, rep- 
resent a major obstacle for the long-time coherent con- 
trol of quantum states. However, in the last years it has 
been realized that dissipation can be extremely useful if 
it can be controlled accurately. Recent experiments have 
shown that strong correlations can be induced by two- 
body losses in ultracold quantum gases ^^. Three-body 
losses can be tailored to generate effective three-body in- 
teractions !3] and to prepare strongly correlated states 
for quantum simulations of color superfluidity j^] , quan- 
tum hall physics [5| or d-wave pairing [6|. Even more, 
dissipation can be used as a universal tool in quantum 
state preparation [7, 8], entanglement generation 9] and 
quantum information processing [10[. These concepts of 
controlling quantum dynamics and transport are partic- 
ularly important for experiments with ultracold atoms in 
optical lattices, where it is possible to address the quan- 
tum system with single-site resolution [ll|, Hjl • An even 
higher spatial resolution has been realized with a focussed 
electron beam, removing atoms one-by-one from the lat- 
tice [ij,[ij] ■ The effects of such a localized particle loss on 
the dynamics of a Bose-Einstein condensate (BEC) have 
been investigated from a nonlinear dynamics viewpoint in 
several papers in the last years, discussing the possibility 
to induce nonlinear structures such as bright breathers 
[isl [l6| , dark solitons [ij] or ratchets [131 ■ These studies 
were based on a mean-field approximation, where the loss 
was introduced heuristically as an imaginary potential. 

In this article we go beyond this approximation and 
investigate the quantum dynamics of ultracold atoms in 
a finite optical lattice with dissipation, which provides a 
distinguished model system for the study of open one- 
dimensional chains. Our analysis is based on a numer- 



ical integration of the full many-body master equation 
and generalized mean-field methods. In section [III we 
present an explicit derivation of the mean-field equations 
of motion, which hold if the many-body state is close to 
a BEC, and generalize this approach to take into account 
higher order correlation functions [1^, [201 ■ I^ particle loss 
is the only source of dissipation, the mean-field equations 
reduce to a non-hermitian Schrodinger equation applied 
previously |16l . [l7[ • While such a non-hermitian descrip- 
tion has been thoroughly studied for single particle quan- 
tum mechanics f2l'|, the applicability to open many-body 
systems is an open issue. 

Two important cases are studied in detail: In section 
mil we analyze how boundary dissipation induces local- 
ization and purifies a BEC. In section IIV[ we consider 
localized loss from a single lattice site, which creates a 
vacancy and leads to a fragmentation of the condensate. 
Remarkably, strong dissipation can suppress the decay of 
the condensate and a coherent dark soliton can be gen- 
erated by properly engineering the dynamics. The tech- 
niques presented here can be directly applied in ongoing 
experiments [Ij, [ij] . 



II. THE MEAN-FIELD LIMIT AND BEYOND 

The coherent dynamics of ultracold atoms in opti- 
cal lattices is described by the celebrated Bose-Hubbard 
Hamiltonian 12211 
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where aj and a'- are the bosonic annihilation and creation 
operators in mode j, J denotes the tunneling matrix el- 
ement between the wells and U the interaction strength. 
We set H = 1, thus measuring energy in frequency units. 
This model assumes that the lattice is sufficiently deep, 
such that the dynamics takes place in the lowest Bloch 
band only. 



In the presence of dissipation, the dynamics is given 
by a master equation in Lindblad form [23i] . 



-i[H,p]+Cp. 



(2) 



Here, we are especially interested in the effects of local- 
ized particle loss, which can be implemented by an elec- 
tron beam [13|, |lj] or by a strongly focussed resonant 
blast laser. Furthermore, phase noise is always present 
in experiments, which degrades the phase coherence be- 
tween adjacent wells and heats the sample [ISjIl^- These 
two processes are described by the Liouvillians [2J, \2uh 
M 
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^phascP ^~2^n'^]P + pf^]- 2nj pfij , (4) 

where 7j denotes the loss rate at site j and k is the 
strength of the phase noise. 

To derive the mean-field approximation, we start from 
the single particle reduced density matrix (SPDM) ajk ~ 
{ajCLk) = tr(a'|afc/5) [1^, |20, [29[ • The equations of motion 
for ajk are obtained from the master equation ([2|). 

^jl^^i^fe = ir {a},jak[H , p] + ia^cikCpj 

= —J {Oj^k+l + 0'j,fe-l — ^j+l,k — (7j-l,k) 
-\-U {(JkkCrjk + ^jkkk — (^jj^jk — ^jjjk) , 

.7j +7fe 
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where we have defined the covariances 
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In the mean-field limit N ^ oo with g — UN fixed, one 
can neglect the variances Ajk^m in Eq- © in order to 
obtain a closed set of evolution equations. This is appro- 
priate for a pure BEC, because the variances scale only 
linearly with the particle number N, while the products 
o'jkO'im scale as N"^. If phase noise can be neglected, i.e. 
K = 0, the equations of motion ^ are equivalent to the 
non-hermitian discrete nonlinear Schrodinger equation 

i-Tfi'k = -J{^^^k+l + VA'-i) + U\i>k\'^^pk - i-^'^l^k (7) 

by the identification cTj^k = ^*j'4'k- This provides a proper 
derivation of the non-hermitian Schrodinger equation, 
which has previously been applied heuristically |15l4l7l |. 
The mean-field approximation assumes a pure BEC 
and is strictly valid only in the limit N ^ oo. To de- 
scribe many-body effects such as quantum correlations 
and the depletion of the condensate for large, but finite 
particle numbers, we gene ralize the Bogoliubov backre- 
action (BBR) method !1£| to the dissipative case, taking 
into account the covariances ^ explicitly. We start with 



the coherent part of the master equation, which yields 
the following evolution equations for the four-point func- 
tions: 
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Again, the interaction hamiltonian leads to higher-order 
correlation functions. To obtain a closed set of evolution 
equations, these function are truncated according to [20| 

(aU„a[, dnfltds) « (djdmd^d„)(d];ds) 



-(d]dmaj;as)(d[a„) + (a{,d„d];ds)(d]d, 
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For a BEC, the six-point function scale as N'^ ^ while the 
error introduced by this approximation increases only lin- 
early with A^. The relative error induced by the trunca- 
tion thus vanishes as 1/iV^ with increasing particle num- 
ber. Close to a pure condensate, the BBR method thus 
provides a better description of the many-body dynam- 
ics than the simple mean-field approximation, since it 
includes the dynamics of higher order methods at least 
approximately. Using this truncation, the coherent part 
of the dynamics is given by 

.d_ 
dt 

^jjn.k+l.n ^j,m.fc — l,n ^j + l^ni.k.n ^j — l,m,/i;,nj 
\U ^^mmknO'jfn ^jjkji^jin f ^jtnnn^kn ^jmkk^kn 
+Ajmfc„ {(Jmm + <^nn — (^kk — f^jj ) J • (10) 

Particle loss and dissipation affect the dynamics of the 
four-point functions as follows 
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FIG. 1: (Color online) Numerical test of the BBR methods 
for a leaky double-well trap with loss in the second well. (a,b) 
Dynamics of the population imbalance (n2 — ni) for two dif- 
ferent values of the loss rate, comparing the BBR approxima- 
tion (solid blue line) to numerically exact results (thick red 
line), (c) Condensate fraction Ao/ntot as a function of time 
and the loss rate 72. (d) Trace distance ([12} between the 
exact rescaled SPDM a(t)/n(t) and the respective BBR ap- 
proximation. In all cases the initial state is assumed to be a 
pure BEC with with equal population and a phase difference 
of TT between the two modes. The remaining parameters are 
J = 10s~\ K = 0,U = 0.5 s"^ and n(0) = 200 atoms. 



The BBR method is especially useful if the many-body 
state is close to, but not exactly equal to a pure BEC. 
In particular, it accurately predicts the onset of the de- 
pletion of the condensate mode. The number of atoms 
in this mode is given by the leading eigenvalue Aq of the 
SPDM aj^k, where the trace of aj^k gives the total num- 
ber of atoms ritot- The ratio Ao/ntot is referred to as the 
condensate fraction [I3, |30| • 

The BBR approach has been extensively tested for 
closed systems in [20] ■ Therefore, we only briefly com- 
ment on the performance of this method in the presence 
of dissipation. Figure [1] shows two examples of the dy- 
namics of a BEC in a leaky double-well trap, comparing 
the BBR approximation (solid blue line) and numerically 
exact results (thick red line). The initial state is assumed 
to be a pure BEC with equal population and a phase dif- 
ference of IT between the two modes. In the case of strong 
dissipation, the BBR approximation predicts the correct 
evolution of the population imbalance (^2 — fii) with an 
astonishing precision. In contrast, significant differences 
are observed for weak losses. This means that the pres- 
ence of particle loss actually improves the performance 
of the BBR method, as the dissipation drives the many- 
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FIG. 2: (Color online) Dynamics of a BEC in a triple- well trap 
with boundary dissipation: (a) atomic density ifijit)), (b) 
total particle number and (c) the condensate fraction Ao/wtot 
for J = 5s-\ 7 = 20s-\ K = 0, [/ = 30 s"^ and n(0) = 60 
atoms, (d) Total particle number after a fixed propagation 
time tfinai = 0.2 s as a function of the interaction strength 
U . Mean-field (— • — ) and BBR ( — ) results are compared to 
numerically exact simulations with a quantum jump method 
averaging over 200 trajectories (thick red line). 



body quantum state towards a pure BEC [3y|. This is 
confirmed by the numerical results for the condensate 
fraction Ao/ntot plotted in Fig. [T](c). A significant de- 
pletion of the condensate is only observed for small values 
of the loss rate 72 . For a further quantitative analysis of 
the accuracy, we compare exact and BBR results for the 
rescaled SPDM <j{t)/ntot{t)- Figurc[l](d) shows the trace 
distance of the exact matrix and the matrix obtained by 
the BBR method. 



d := 2*'^(I'^bbr/"bbr - CTox/^cxI), 



(12) 



as a function of time for different values of 72 . For suf- 
ficiently large dissipation, one observes that the distance 
approximately vanishes for all times. In this regime the 
quantum dynamics is faithfully reproduced by the BBR 
approximation. 



III. BOUNDARY DISSIPATION 

We first analyze the effects of boundary dissipation 
with a focus on small systems for which numerically ex- 
act solutions of the many-particle dynamics are still pos- 
sible, for instance by the quantum jump method |23ll3ll|. 
A comparison to numerically exact results for these ex- 
amples provides another test of performance of the BBR 
approach. 



We consider the decay of an initially pure, homoge- 
neous BEC in a triple-well trap with boundary dissipa- 
tion. Figure [2] (a) and (b) show the evolution of the 
atomic density and the total particle number for strong 
inter-atomic interactions U — 30 s^^. One observes a 
fast decay of the atoms at the outer sites while the pop- 
ulation at the central site is remarkably stable. This is 
confirmed by the evolution of the total particle number, 
which rapidly drops to about one third of its initial value, 
where it saturates for a long time. This is a consequence 
of the dynamical formation of a discrete breather at the 
central site, which is an important generic feature of non- 
linear lattices. Generally, discrete breathers, also called 
discrete solitons, are spatially localized, time-periodic, 
stable excitations in perfectly periodic discrete systems 
|33 - l35| . They arise intrinsically from the combination of 
nonlinearity and the discreteness of the system. In the 
presence of boundary dissipation, these excitations be- 
come attractively stable such that the quantum state of 
the atoms will converge to a pure BEC with a breather- 
like density for a wide class of initial states. Once a 
discrete breather is formed, it remains stable also if the 
dissipation is switched off. The crucial role of strong in- 
teractions is illustrated in Fig. [2](d), where the residual 
atom number after ifinai — 0.2 s of propagation is plotted 
as a function of the interaction strength. The particle 
number increases for large values of U to ntot(^finai) ~ 20 
due to the breather formation. 

For strong interactions a simple mean-field approxima- 
tion fails. It strongly underestimates the residual par- 
ticle number as it predicts that discrete breather are 
formed only for stronger losses. In contrast, the BBR 
results agree well with the many-particle simulation even 
for large values of U. We thus conclude that quantum 
fluctuations facilitate the formation of repulsively bound 
structures. Furthermore, a mean-field approach cannot 
account for genuine many-body features of the dynam- 
ics. Figure [5] (c) shows the evolution of the condensate 
fraction Ao/ntot, where Aq is the leading eigenvalue of 
the SPDM [3C|. In the beginning, interactions lead to a 
rapid depletion of the condensate. On a longer time scale, 
however, dissipation restores the coherence and drives the 
atoms to a pure BEC localized at the central lattice site 
|3a | . The BBR approach faithfully reproduces the deple- 
tion and re-purification but additionally predicts unphys- 
ical temporal revivals. This example thus demonstrates 
the strength but also the limitations of this method. 

The decay dynamics of the discrete breather state is 
further analyzed in Fig.|31 The total atom number ntot(^) 
decreases rapidly until the discrete breather is formed at 
t « 0.2s. Afterwards the decay is much slower and clearly 
non-exponential. In both regimes, one can calculate the 
evolution of ritot (t) approximately, starting from the re- 
lation ntot = —"f{ni + n^). Initially, all sites are filled 
homogeneously, ni —713= ritot/S, such that the total 
particle number decays as 




FIG. 3: (Color online) Decay of a discrete breather state for 
J = 5s"^ 7 = 20s~\ «: = and U = lOs"^ Numerical 
results calculated with the BBR method ( — ) are compared 
to the analytic estimates (|13p and (|14l) . respectively ( ). 



When the discrete breather is formed, the population of 
the outer wells is given by ni = ria = J^/{U^ntot) in 
first order perturbation theory. The atom number then 
decays as 



ntot(i) « sJnl^-4jJH/U^ 



(14) 



ntot(i)«ntot(0)e-27/3*. 



(13) 



where n^h is the number of atoms bound in the discrete 
breather state. Both approximations are compared to the 
BBR simulation results in Fig. [31 assuming a breather 
with ridb = 13 atoms. One observes an excellent agree- 
ment in the both regimes, i.e. an exponential decay for 
very short times {t < O.l s) and an algebraic decay when 
the discrete breather is formed. The transition between 
the linear and nonlinear decay takes place at t « 0.2 s. A 
deviation from the algebraic decay (fT4| for the discrete 
breather is observed only for very long times when the 
atom number is very small such that the simple pertur- 
bative estimate for rii_3 is no longer valid. 



IV. LOCALIZED LOSS 

Recent experiments with ultracold atoms have demon- 
strated an enormous progress in spatial addressability us- 
ing specialized optical imaging systems [11,, 12] or a fo- 
cussed electron beam [lj,[l3|- Especially the latter exper- 
iment allows to manipulate a Bose-Einstein condensate in 
an optical lattice dissipatively with single-site resolution. 
In the following, we study the quantum dynamics in a 
finite lattice of 11 sites with closed boundary conditions 
and loss occurring from the central site only, which leads 
to remarkably different decay as in the case of boundary 
dissipation studied above. 

A remarkable feature of the quantum dynamics is il- 
lustrated in Figure 21 showing the results of a BBR sim- 
ulation for an initially pure homogeneous BEC. For a 
modest loss rate 7 = 20 s~^, atoms tunnel to the central 
site where they are dissipated with a rate 7, such that 
the BEC decays almost homogeneously. On the contrary, 
stronger losses (7 — 100 s~^) lead to a formation of a sta- 
ble vacancy. The central site is rapidly depleted, but the 
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FIG. 4: (Color online) Formation of vacancies by localized 
loss at the central lattice site. (a,b) Evolution of the atomic 
density {fijit)) (colorscale as in Fig. [2|. (c,d) Final value of 
the total particle number after a fixed propagation time iflnai 
as a function of the loss rate 7, without (solid lines) and with 
strong phase noise (dashed lines, k — 50 s~^). The remaining 
parameters are J — 5s~^, U = Q.2s~^ and n(Q) = 1000 
particles (a-c) and U = 2s~^ and n(0) = 50 particles (d). 
The dynamics has been simulated with the BBR (thin blue 
lines) and the quantum jump method (thick red lines). 



atoms in the remaining wells are mostly unaffected. Thus 
one faces the paradoxical situation that an increase of the 
loss rate can suppress the decay of the BEC. Two effects 
contribute to this counterintuitive behavior: (i) The ab- 
sorbing potential suppresses tunneling to the leaky lattice 
site. This effect is present also in the linear case and can 
be explained by an analogy to wave optics [l|: A large 
mismatch of the index of refraction leads to an almost 
complete reflection of a wave from a surface. This is true 
for an imaginary index describing an absorption as well 
as for a real index, (ii) A dark breather stabilizes the va- 
cancy and prevents the flow of atoms to the central site. 
This nonlinear structure remains stable also if the dissi- 
pation is reduced or switched off afterwards (cf. [33-35| 
for a discussion of the stability of breathers). 

The suppressed decay of the BEC is further illustrated 
in Fig. m (c,d), where the residual atom number after a 
fixed propagation time is plotted as a function of the loss 
rate 7. The coherent output of the system, i.e. the num- 
ber of lost atoms, assumes a maximum for a finite loss 
rate 7c. This maximum is reminiscent of the quantum 
stochastic resonance discussed in [36|. In the following 
we will estimate the value of 7c by determining a lower 
bound for 7 for the dynamical breather formation. As 
a single (both bright and dark) breather exhibits a pro- 
nounced population imbalance between the central site 
and the neighboring sites, we estimate 7c by matching 
the timescales of dissipation tjj = 2/7 and tunneling r,/, 
i.e., T£) = Tj. For smaller values of 7, atoms can tun- 
nel away from the leaky lattice site again before they are 
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FIG. 5: (Color online) Coherence of a vacancy generated by 
loss from the central site: (a,b) atomic density, (c,d) scaled 
eigenvalues Am/ntot of the SPDM and (e,f) phase coherence 
g^^' between the two BEC fragments. Parameters are the 
same as in Fig. 0] with 7 = 100 s~^ and «; = 0. Results of a 
quantum jump simulation are plotted as thick red lines, BBR 
results as thin blue lines. 



lost, while for larger values of 7 a population imbalance 
can form. From Eq. ([7]) we read tj = 1/(2J) where the 
factor 1/2 accounts for atoms tunneling from two sites to 
the leaky site. Hence, the critical loss rate is estimated as 
7c = 4 J. We find good agreement of our qualitative esti- 
mate for 7c (dotted vertical lines in Fig.|3](c,d)) with the 
dip in the total particle number. An important quantity 
for the breather formation and stability is the effective 
nonlinearity of the system A = Untot{t)/2J, which, due 
to particle loss, is time-dependent. Strikingly, though A 
depends on the interaction strength U (which is differ- 
ent in Fig. |l](c) and (d)), the fairly good estimate 7c is 
independent of U. 

Figure |4] (d) shows the respective results for a triple- 
well trap with loss from the central site. A comparison 
of the BBR approximation to a numerically exact many- 
particle simulation shows a good agreement for all values 
of 7. Phase noise suppresses decay as it effectively de- 
couples the lattice sites. Thus, only the atoms initially 
loaded at the leaky lattice site decay as e"'''*, while the 
other atoms remain at their initial positions. With in- 
creasing loss rate 7, the number of atoms lost from the 
trap approaches ~ n{0)/M as shown in Fig. |3](c,d). 

The previous reasoning suggests to use dissipation as a 
tool to coherently engineer the quantum state of a BEC 
in an optical lattice. Mean-field theory predicts that dis- 
sipation can be used to efficiently create a coherent dark 
soliton 17| , but cannot assert the coherence of the final 
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FIG. 6: (Color online) Generation of dark solitons using loss 
imprinting at a rate 7 = 100 s~^ at the central site (a,c) and 
phase imprinting in the lower half of the lattice (b,c), both 
for times t < 0.1s only. Shown are the atomic density (left, 
colorscale as in Fig.[2| and the scaled eigenvalues of the SPDM 
(right) calculated with the BBR method. Parameters are J = 
5s"\ U = O.ls"^ K = and n(0) = 1000 particles. 



state as discussed above. The results of a BBR and a 
quantum jump simulation of the many-body dynamics 
shown in Fig. [S] reveal the limitations of the phase coher- 
ence of a soliton generated by local dissipation. The up- 
per panels (a^b) show the rescaled eigenvalues Am/n-tot of 
the SPDM [30|. One observes that there are two macro- 
scopic eigenvalues approaching 1/2, while all remaining 
eigenvalues vanish approximately. This proves that the 
dissipation generates a fragmented BEG consisting of two 
incoherent parts rather than a single BEG with a soli- 
tonic wavefunction. The BBR simulations correctly de- 
scribe the fragmentation of the condensate, but predict 
temporal revivals of the coherence which must be consid- 
ered as artifacts of the approximation. Experimentally, 
one can test the coherence by the interference of the two 
fragments in a time-of-fiight measurement. Figure[5](e,f) 
shows the coherence 
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between the wells £ and m. One clearly observes the 
breakdown of phase coherence between the two conden- 
sate fragments. 

In order to overcome the loss of coherence, one can, 
however, engineer the many-body dynamics. Figure [6] il- 
lustrates the generation of dark solitons comparing three 
different strategies. If the dissipation is switched off after 
the generation of a vacancy at t = 0.1s, the condensate 
remains pure for long times. However, the vacancy is not 
stable but decays into two dark solitons traveling out- 



wards [iTf, where they are reflected at the boundaries. 
The effects of a phase imprinting, which is an established 
experimental method [33|, are shown in Fig. |6] (b). A 
local potential is applied to the lower half of the lattice 
for t < 0.1s imprinting a phase difference of it. Again 
coherence is preserved but the generated solitons travel 
outwards. A coherent and stable dark soliton can be en- 
gineered by combining both methods, as shown in Fig. [S] 
(c). The generated dark soliton stays at its initial posi- 
tion and remains coherent over a long time. 



V. CONCLUSION 

We have discussed the influence of localized particle 
dissipation on the dynamics of a finite one-dimensional 
Bose-Hubbard chain, which describes a Bose-Einstein 
condensate in a deep optical lattice [Ij, [lj|. Starting 
from the many-body master equation, we have derived 
the mean-field approximation and the dissipative Bogoli- 
ubov backreaction method, which allows a consistent cal- 
culation of the depletion of the condensate. 

Two important special cases have been studied in de- 
tail. Particle loss at the boundary leads to localization 
and the formation of coherent discrete breathers. Sur- 
prisingly, dissipation together with interactions can re- 
purify a BEG. A striking effect of localized loss is that 
strong dissipation can effectively suppress decay and in- 
duce stable vacancies. The decay shows a pronounced 
maximum for intermediate values of the loss rate, when 
the timescales of the dissipation and the tunneling are 
matched. Gombined with an external potential, these 
effects can be used to generate stable coherent dark soli- 
tons. These examples show that engineering the dissipa- 
tion is a promising approach for controlling the dynamics 
in complex quantum many-body systems. 

Ultracold atoms provide a distinguished model system 
for the dynamics of interacting quantum systems, such 
that the effects discussed in the present paper may be 
observed in different systems, too. In particular, quan- 
tum transport of single excitations driven by local dissi- 
pation has recently been studied in a variety of physical 
systems ranging from spin chains '38* to light-harvesting 
biomolecules [32,]. On the other hand, it has also been 
shown on the mean-field level that nonlinear excitations 
such as discrete breathers play an important role for 
quantum transport in these systems (cf. [33, 3J, |40| and 
references therein) . Thus it is of general interest to fur- 
ther explore the regime which interpolates between the 
nonlinear mean-field dynamics and the many-body quan- 
tum dynamics in the spirit of the work presented here. 
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